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Abstract 

One key feature of isogeometric analysis is that it allows smooth shape functions. Indeed, when 
isogeometric spaces are constructed from p-degree splines (and extensions, such as NURBS), they 
enjoy up to continuity within each patch. However, global continuity beyond on so-called 

multi-patch geometries poses some significant difficulties. In this work, we consider planar multi¬ 
patch domains that have a parametrization which is only at the patch interface. On such 
domains we study the /i-refinement of -continuous isogeometric spaces. These spaces in general 
do not have optimal approximation properties. The reason is that the C^-continuity condition 
easily over-constrains the solution which is, in the worst cases, fully locked to linears at the patch 
interface. However, recently m has given numerical evidence that optimal convergence occurs for 
bilinear two-patch geometries and cubic (or higher degree) splines. This is the starting point 

of our study. We introduce the class of analysis-suitable geometry parametrizations, which 

includes piecewise bilinear parametrizations. We then analyze the structure of isogeometric 
spaces over analysis-suitable parametrizations and, by theoretical results and numerical testing, 
discuss their approximation properties. We also consider examples of geometry parametrizations 
that are not analysis-suitable, showing that in this case optimal convergence of G^ isogeometric 
spaces is prevented. 


1. Introduction 

Thanks to the use of smooth B-splines and NURBS, isogeometric methods [Eiini have revi¬ 
talized the interest for the use of smooth approximating functions for the numerical solution of 
partial differential equations. Advantages with respect to G^ finite element methods are improved 
accuracy and spectral properties giiiaiisiiii] and the possibility to directly discretize differential 
operators of order higher than 2. In the literature there are indeed many examples of isogeometric 
methods for 4^^ order differential problems of relevant interest, such as Kirchhoff-Love plates/shells 
[3 [71 [23], the Cahn-Hilliard equation M, and the Navier-Stokes-Korteweg equation m- 

Since higher dimensional spline spaces possess a tensor-product structure, the representation of 
domains that have a complex geometry is non-trivial. In this paper we focus on multi-patch repre¬ 
sentations. While the implementation of C^-continuity over multi-patch domains is well understood 
(see e.g. for strong and [5] for weak imposition of the conditions), C^-continuity is 
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not. Several studies have tackled the problem of constructing function spaces of or higher order 
smoothness. A first attempt to compare different ways to impose -continuity in isogeometric 
analysis was presented in [29]. We also refer to [inmaiisi for smooth constructions for spline 
spaces and 12 ^ 1110 ] for triangulations, which can be seen as an alternative to the classical B-spline 
based isogeometric framework. Nevertheless, the construction of smooth isogeometric spaces with 
optimal approximation properties on complex geometries is still an open and challenging prob¬ 
lem. This is related to the problem of finding parametrizations of smooth surfaces having complex 
topology, which is a fundamental area of research in the community of Computer Aided Geometrid 
Design (CAGD) over the last decades. 




Figure 1: Two possible parametrization schemes: away from extraordinary points (left) and C° everywhere (right). 

We review two different strategies for constructing smooth multi-patch geometries and corre¬ 
sponding isogeometric spaces. One possibility is to adopt a geometry parametrization which is 
globally smooth almost everywhere, with the exception of a neighborhood of the extraordinary 
points (or edges in 3D), see Figure(left). The other possibility is to use geometry parametriza¬ 
tions that are only at patch interfaces, see Figure(right). The first option includes subdivision 
surfaces im and the T-spline construction in [3H] and, while possessing attractive features, they 
seem to possess optimal approximation properties for some configurations, see m, but in general 
lack accuracy HSIES]. In our work we consider the second possibility, corresponding to the right 
part of Figure 

The construction of isogeometric functions over a parametrization can be interpreted con¬ 
veniently as geometric continuity of the graph parametrization. Bilinear multi-patch parametriza¬ 
tions of a planar domain have been analyzed in 13 EH, where it is shown that there exists a minimal 
determining set with local degrees of freedom for the space of (mapped) piecewise polynomial func¬ 
tions, with global continuity, if the polynomial degree is high enough (4 if some additional 
conditions are fulfilled, 5 in general). The recent preprint [28| generalizes the previous results, by 
using advanced homology techniques, to arbitrary parametrizations, allowing both triangular and 
quadrilateral patches. 

In the work m , the authors consider splines instead of polynomials within each patch, construct 
a basis, analyze the space dimensionality of some configurations and, for the first time, perform 
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numerical tests to evaluate the order of convergence when each patch is /i-refined. We recall that 
within the isogeometric framework, the concept of /i-refinement, equivalent to knot insertion, is 
one of the three constructions to increase accuracy of the spline spaces (see [ED- Remarkably, m 
gives numerical evidence of optimal convergence for splines of degree 3 (or higher), on a two- 
patch bilinear geometry. The authors also show an example of over-constrained isogeometric 

spaces, corresponding to a two-patch non-bilinear geometry parametrization. We refer to the latter 
situation as locking. Our work develops the underlying theory. While the previous papers give 
explicit charaterizations in the form of minimal determining sets mEzi or basis constructions m 
EH, we use an implicit characterization of the continuity conditions and derive from it information 
on the structure of the isogeometric space. As in |21| . our interest is in the impact of /i-refinement. 
We study /i-refinement for arbitrary degree and regularity, both theoretically and numerically, and 
point out its performance depending on the geometry parametrization. 

We set up our notation in Section]^ In Section]^ we define the class of analysis-suitable (AS) 
geometry parametrizations, which includes the bilinear ones and the extensions presented in 
Section 3.4 of EH- Then, in Section |5.l| we study the structure of isogeometric spaces over 


AS two-patch geometries. Here, we give an explanation of the optimal convergence of p-degree 
isogeometric functions, having G^ continuity across the patch interface and up to continuity 

within the patches. Furthermore, we discuss why G^ locking occurs for G^~^ continuity within 
the patches. Note that in this paper we do not derive explicit error estimates, which will be 
the topic of a future paper. In Section 5^ we analyze G^ isogeometric spaces constructed over 
more general geometry parametrizations and conclude that /i-convergence is suboptimal beyond 
AS G^ geometries. The extensions to surface domains and to NURBS are briefly discussed in 
Sections and respectively. Numerical tests on two- and multi-patch domains are reported in 
Section There we present a significant example: the multi-patch parametrization of a smooth 
simply-connected planar domain. The question of existence and construction methods for AS G^ 
multi-patch parametrizations of arbitrary geometries remains to be studied. We summarize our 
results and draw conclusions in Section |9l 


2. Planar multi-patch spline parametrizations and isogeometric spaces 

Given an interval or a rectangle cj, we denote by Sr(cj) the spline space of degree p (in each 
direction) and continuity G^ at all interior knots. The knot mesh in the parametric domain is 
assumed to be uniform, with mesh-size h (which is not explicitly indicated in the notation) and 
interior knot multiplicity p — r. We write Sr instead of Sr(oj) when the domain cv is obvious from the 
context. We allow r > p, which stands for G^ continuity, that is, the case of global tensor-product 
polynomials on cj. In this case we also use the notation = Sp{uj) = <s^+iM =.... 

Consider a planar multi-patch domain of interest 

U ... U C (1) 

where the closed sets form a regular partition (with disjoint interior). For simplicity, we do 
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(2) 


not allow hanging nodes. Each 


is assumed to be a spline patch, that is 


F® : [0,1] X [0,1] = ^ fl®. 


where G Sr{^) x 5r (f^). We assume 

r > 1 (3) 

{r > p means we have Bezier patches and we assume the parametrizations are not singular, 
i.e., for all i and for all (u^v) G fl. 


det 


DuF^^{u,v) DyF^^\u,v) 


^0. 


( 4 ) 


For the sake of simplicity we do not consider more general configurations, e.g., non-uniform knot 
meshes, different degree or continuity parameters in each parametric direction, different continuity 
at different knots, or different spline spaces for different patches. Indeed, our simple configuration 
already presents the key features and difficulties we are interested in. 

We assume global continuity of the patch parametrizations. This means the following. Let us 
fix r = n when using the superscript as (i, j) in the paper, we assume implicitly 

that i and j are such that is not a point or an empty set. Let F^^), F^^) be given such that 

F(^) : [-1, 0] X [0,1] = = fl®, 

: [0,1] X [0,1] = ^ 


where (F(-^))-1 oFW and are linear transformations (combinations of a translation, 

rotation and symmetry). Moreover, the parametrizations agree at ix = 0, i.e., there is a Fq : [0,1] ^ 
with 

r = {Fo(n) = F(^)(0,n) = F(^)(0,n),nG [0,1]}. (6) 

An example is depicted in Figure 



Figure 2: Example of the general setting of 
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Remark 1. The domain Q = U ... U can be endowed with a spline manifold structure as 
defined, e.g., in naiiaEf. In the framework of m, each pair of adjacent subsets is 

naturally associated with a chart [—1,1] x [0,1] = U through the maps and 


Definition 1 (Isogeometric spaces). The isogeometric space corresponding to Sr and Tt is given as 


v = {, 

^ : n ^ R such that 0 o G 5^(12), i = 1,..., Al} • 

(7) 

Furthermore we have 


yo = vnc°(fi), 

(8) 

and 


yi = yncHu). 

(9) 


The graph D C X R of an isogeometric function 0 : f2 —)■ R is naturally split into patches S* 
having the parametrizations 


ph) 

n(i) 


: [0,1] X [0,1] = fl ^ E® 


( 10 ) 


where = cf) o 

In order to analyze the smoothness of an isogeometric function along one interface F = = 

n we introduce 


f(-^) 

j,(i) 

f(-^) 

giR) 


: [-1,0] X [0,1] = sW = 


: [0,1] X [0,1] = 


( 11 ) 


where are dehned obviously as extensions of ([^, see Figure]^ Continuity of cj) is implied 

by the continuity of the graph parametrization, which we assume and set to 


9 o{v) = g^^\0,v) = g^^\o,v), 


( 12 ) 


for all V G [0,1], analogous to (§. 


3. isogeometric spaces 

If an isogeometric function is within each patch (condition Q), with non-singular parametriza¬ 
tion (condition Q) and is globally continuous (condition ([T^), then it is globally if and only 
if there exists a well defined tangent plane at each point of the interfaces H Focusing on 

one interface, with notation 0, the tangent planes from the left and right sides are formed by 
the two pairs of vectors 


I1„fW(0,u) 

and 

-D^Fo(u) 

as well as 

F>„F(^)(0,u) 

and 

Zl^Fo(u) 

Dug^^\^.v) 


Dvgo{v) 


Dug^^\Q,v) 


Dvgo{v) 
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Figure 3: Example of the general setting of (11). 


respectively. These are three different vectors (the vector tangent to n is in common) that 
form a unique tangent plane, i.e. they are coplanar, if and only if they are linearly dependent. In 
other words the isogeometric function 0 is on U if and only if 


n„F(^)(0,^;) n„F(^)(0,^;) n,Fo(^;) 

Dug^^Ko,v) Dug^^\0,v) Dygo{v) 


for all V G [0,1]. In the context of isogeometric methods, the domain Q and its parametrization 
are given at the first stage, then ( [T^ is the condition on the isogeometric function in parametric 
coordinates (i.e., and that gives continuity of the isogeometric function in physical 

coordinates. 


In CAGD literature, condition (13) is named geometric continuity of order 1, in short and 
is commonly stated as in the following Definition (see, e.g., [a ESI IS]). 


Definition 2 (G^-continuity at H Given the parametrizations g^^\ g^^^ as 

in (U), ([n]), fulfilling Q , Q and ( [T^ , we say that the graph parametrization is at the interface 
n if there exist : [0,1] M, : [0,1] ^ M and /3 : [0,1] —>■ M such that for all 
V e [0,1], 

a^^\v)a^^\v) > 0 (14) 


and 


a^^\v) 


DuF(^\0,v) 

Dug^^Ho,v) 


— a^^\v) 


Zl„F(«)(0,n) 

Dug^^\^,v) 


+ ^(n) 


L>„Fo(n) 

Dvgo{v) 


= 0 . 


(15) 


The sign condition on and forbids, on general surfaces, the presence of cusps. 
However, for the graph of a function it is obvious. 

For our study it is useful to express G^ continuity as in Definition since the coefficients 
and /3 play an important role. Since the first two equations of © are linearly independent. 
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Q/(^) and (3 are uniquely determined, up to a common multiplicative factor, by and 
i.e. from the equation 


a^^\v)DuF^^\0,v) - a^^\v)DuF^^\0,v) + p{v)DyFo{v) = 0 . 

Precisely, we have the following proposition which can also be found in [SaESlES]. 


(16) 


Proposition 1. Given any in the setting of Section^ then (16) holds if and only if 

= 7(e a for S G {L,R}, and (3{v) = '-y{v)f3{v), where 


a^^\v) = det [ DuF^^1{0,v) DyFo{v) 


(17) 


/3(v) = det DuF^^\0,v) DuF^^\0,v) 


(18) 


and 7 : [0,1] ^ M is any scalar function. In addition, ^{y) ^ 0 if and only if (14) holds. Moreover, 
there exist functions for S G {L,R}, such that 


P{v) — 


(19) 


Proof. Obviously (@ determines and ^{v) up to a common factor 7 ('L’) and, since 

the two vectors DuF^^\0,v) and DyFo{v) are linearly independent because of (W, a^^\v) 


and I3{v) are uniquely determined up to a common factor 7 ('L’) by (16). We have 


det 


I1„F(^)(0,^;) I1„F(^)(0,^;) Il.Fo(^;) 

Z1„F(^)(0, v) ■ Gk 71„F(^)(0, v) ■ Bk DyFo{v) ■ 


for all = 1, 2, ek being the canonical base vectors in I 
above determinant (along the third row) we end up with 


= 0 , ( 20 ) 

Using the Laplace expansion of the 


a^^'>iv)DuF^^\0,v) - d^\v)DuF^^'>{0,v) + P{v)D^Fo{v) = 0, 

where and j3 are as in 0 and ( [T^ . In a similar way ( [T^ yields 

a^^\v)Dug^^\0,v) - a^^\v)Dy^\0,v) + j3{v)D^gQ{v) = 0 . 

Setting 


we have 


u{v) = 


det 


DyFo^v) • 62 
—DyFo{v) • 61 


t{v) = 


DvFq{v) 

|^.Fo(n)||2’ 


iy{v) t(v) 


= 1 and iy(v) • t(v) = 0. 


The existence of S G such that (19) holds is obvious, because of (14). Obviously, 

for S G are not unique. A specific choice, which is of interest, is the following 


(S). . ^ Il„F(^)(0,n)-Zl,Fo(n) 
^ ^ ^ P.Fo(n)P 


( 21 ) 
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Indeed, using the expansion 


\/v G 


V = det 


V t{v) 


iy{v) — det 


V u{v) 


t{v), 


and then 


det 


V 


{L) y{R) 


= det 


V 


{L) 


u{v) 


det 


V 


(R) 


t{v) 


det 


t{v) 


det 


u{y) 


gives (19) and (21) by choosing ?;), for S G 


□ 


Remark 2. If (15) holds, then there exist eoeffieients G 1]), for S G {L,R}, and 

(3 G 5r^([0,1]). Indeed, this follows from Proposition^ seleeting 7 = 1. See also m- 

Summarizing, in the context of isogeometric methods we consider 12 and its parametrization 
given. Then for each interface and /3 are determined from (16) as stated in Proposition 

It should be observed that for planar domains, there always exist and /3 fulfilling (16). 

This is not the case for surfaces, see Section Then, the continuity of isogeometric functions 
is equivalent to the last equation in , that is 

a^^\v)Dug^^\(),v) - a^^\v)Dug^^\^.,v) + l5{v)D^go{v) = 0 ( 22 ) 


for all V G [0,1]. 

We end this section by a statement of the equivalence between continuity of the isogeometric 
function and continuity of its graph parametrization. It is formalized and presented in its most 
general form for arbitrary continuity and dimension in d- The use of continuous functions 
over unstructured mesh partitions is well known in the isogeometric community, see e.g. iiEiiEa 
EHEZIEHI. We give a detailed proof of the statement here, in the framework of Proposition since 
this will serve for the next steps of Section 

Proposition 2. An isogeometrie funetion 0 G V belongs to if and only if its graph S is G^ 
eontinuous on eaeh interfaee H 


Proof Consider a graph interface n and the corresponding P = P^^d) = p qU) ^ Let 
for S G {L,R} such that (16) holds. Define the vector on P such that 


d^^KFoiv) = \ DuF^^\0,v) D^Fo{i 




a 




(23) 


The vector is transversal to P, i.e. linear independent to DyFo{v), since 

1 


det 


d(^)oFo(^;) D,Fo{i 


a('^)(v) 


det 


Z)„F(^)(0,^;) n,Fo(^ 


l(v] 


^0. 


We have = d. Indeed, by using © and we get 

o Fo('r’) — o Fo{v) 

= a^^\v)DuF^^\0,v) - a^^\v)DuF^^\0,v) + - a^^\v)^^^\v))D^Fo{v) 


= 0 . 

































Futhrermore, d is not in general unitary and it is continuous, since Fq is at least from (§. For 
S G let (j) = o on then the derivative from side S in direction d fulfills 


\/^^Uix,y)-d{x,y) = [Dug^^\o,v) D^goiv)] \ DuF(^\0,v) D^Fo{v 


-1 


doFo(^;), 


(24) 


for all {x,y) G F, where here and in what follows we implicitly assume the relation Fo(n) = {x,y). 
We obtain directly from the definition of that 


DuF(^\0,v)D^Fo{i 


1-1 ^ ^ 1 

1 

■doFo{v) = 

-1 



(25) 


Substituting (25) back to (24) we then obtain 


,). d(..,) = (2f>) 

for S G {i?, L}. Therefore ' d{x^y) — ' d{x^y) if and only if (by ([^) 

Du9^^\Q,v) - P^^\v)Dvg(){v) ^ Dug^^\Q,v) - P^^\v)Dvgo{v) 

That, after multiplying both sides by a^^\v)a^^\v) and using ©, is equivalent to ( [^ . □ 

Remark 3. ^5 a consequence of Proposition\^ isogeometric spaces over planar multi-patch 
parametrizations fit in the isoparametric framework. 

4. Analysis-suitable parametrizations 

At each interface F = n given any regular and orientation preserving 

as in ([^, there exist coefficients and /3, with a^^\v)a^^\v) > 0, such that (16) 

holds. Then (22) expresses continuity of isogeometric functions in terms of and /3. 


Optimal approximation properties of the isogeometric space on hold under restrictions on 
and /3, i.e. on the geometry parametrization. This leads to the definition below. 

Definition 3 (Analysis-suitable G^-continuity). F^^^ and F^^) are analysis-suitable -continuous 
at the interface F (in short, AS G^(F) or AS G^) if there exist G 7^^([0,1]) 


such that (16) and (19) hold, that is for all v G [0,1] 

a^^\v)DuF^^\0^v) — a^^\v)DuF^^\o^v) = — a^^\v)f3^^\v))DyFo{v). 

The degrees of the functions and (3 were also studied in the context of G^ interpolation 

of a mesh of curves in [33[[3lj. There, the same degrees where derived for interpolations using cubic 
patches. 

Remark 4. As in Remark we observe that G^ isogeometric spaces over AS G^ multi-patch 
parametrizations fit in the isoparametric framework. 
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The class of AS parametrizations contains the bilinear ones and more 


Proposition 3. Any G x and G x are AS G^- 

continuous atT. Moreover, for any G ^^^([O,!]) strietly positive G ^^([0,1]) 

there exist F(^) G and F(^) G fulfilling 


Proof The statement follows directly from Proposition and (21). 


□ 


Remark 5. The elass of AS parametrizations is wider than only bilinear. We will show some 
examples later in Section where F^^^ and F^^) are higher order polynomials at {0} x [0,1]. 


5. Two-patch geometry 


In this section we analyze the two-patch geometry. This is the simplest geometric configuration 
that allows us to focus on the C^-continuity constraints that are associated with each patch inter¬ 
face. For the space of isogeometric G^ functions, that is defined in 0-®, optimal convergence 
under /i-refinement is known since the results in apply directly. (7^-continuity constrains 

traces and transversal derivatives at F of functions 0 G V^. Therefore the approximation properties 
of the space follow from the ones of traces of functions and transversal derivatives of functions 
of at F. Let for S G {L,R} such that (16) holds. We define the transversal 

vector d to F as in ( [23| ) and introduce the space of traces and transversal directional derivatives on 
F 

= {r 3 {x,y) ^ [4>{x,y),V4>{x,y) ■ d{x,y)], such that ^ G } (27) 

and its pullback 

Vp = {[0, V0 • d] o Fo, such that 0 G V^} . (28) 

With this choice we have that [0, V(/) • d] o Fq = [go^9i] is equivalent to 


g^^Hu,v) = goiv) + {P^^'>{v)go{v) + {v)gi{v)) u + 0{u^ 


(29) 


thanks to (26). 


Remark 6. The transversal veetor d depends on so do Vp and Vp. The 

funetion traee in Vp fulfills 0 o Fq G Sr, this is not true for the transversal derivative Vfi ’ d whieh 
in general is a rational funetion of some degree higher than p (a similar situation is analyzed in 

m)- 

5.1. AS G^ -eontinuous two-pateh geometry 

The next result gives the key properties of the space Vp (defined in (28)) in the case of AS G^ 
parametrizations. 

Theorem 1. Let Pt = U and let F^^^ and F^^^ be AS G^ at the interfaee F = U 

Then 

5f+i([0,l])x5r'([0,l])c)?. (30) 
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Proof. For linear and (3^^^ {S G {L, i?}) and [go^gi] G 1]) x Sr ^([0,1]), we have 


g(^\u,v) = go{v) + {^^^\v)go{v) P a^^\v)gi{v)) u G 


Then the statement is a direct consequence of (29). 


(31) 

□ 


Theorem guarantees that the trace space for the function value {0 o Fq : 0 G V^} includes all 
splines of degree p and regularity at least r + 1, and independently the trace space for the transversal 
derivatives of the function {{\/<f-d)oFo : 0 G V^} includes all splines of degree p—1 and regularity 
at least r. This suggests that Vp enjoys optimal approximation order, and consequently for the 
whole space when r < p — 1. 

However, if r = p — 1 and the parametrization is not trivial, the space Vp suffers of locking, 
that is, /i-refinement does not improve the approximation properties for 0 and V(j)-d independently. 
The following theorem gives some understanding of this phenomenon. 


Theorem 2. Let Pi — U F^^) he AS at the interfaee F = H and 

r = p — 1. Furthermore, let ctnd Gi be two spaees sueh that 

GoxGiC (32) 


If either (3^^^ ^ 0 or ^ 0 then the dimension of Go is independent of h. If and are 
linearly independent, then the dimension of Gi is independent of h. 


Proof Let [po^O] = [0, V(/) • d] o Fq and assume that is a linear function not identically zero. 
Furthermore assume that if there exists a, vq ^ [O?!]? such that l3^^\vo) = 0, then from the 
assumption ^ 0. Using (29), i.e.. 


g^^Hu,v) = go{v) + I3^^\v)go{v)u + O{u^ 


and using g^^\u, v) G then po is a spline in 5^_]^([0,1]) and also pg is therefore go 

is in fact a degree p global polynomial. The same conclusion holds when there exists a 'L’q ^ [0,1] 
such that P^^\vo) = P^^\vo) = 0 but vo is not a knot of the spline space 5^_]^([0,1]) under 
consideration. If, instead, — 0 and 'L’q is a knot of the spline space Sp_i{[0, 1]), 

then po has in general p—1 continuous derivatives at vo- In conclusion if po G Go then it has to be 
piecewise polynomial of degree p on at most two intervals, and is in C^“^([0,1]) globally. Hence 
the dimension of Go is at most p + 2. 

Take now [0,pi] = [0, V(/) • d] o Fg and assume and are linearly independent, that is. 


their linear combination gives the whole space ^^^([0,1]). Now (29), i.e.. 


g^^\u, v) = a^^^{v)gi{v) u + 0{u^), 


AS) 


after taking the derivative in the variable u and evaluating at = 0 gives 
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and, by suitable linear combination for S = L and S = R yields 


Since g^^\u,v) G gi G <Sp_]^([0,1]). Similarly 

C'L^5^^^(0,e + = vgi{v), 

which means vgi G 5^_x([0,1]). in conclusion gi G tS^Zi ([0,1]), therefore the dimension of Qi is at 
most p. □ 

When = (3^^^ — 0 and for all v G [0,1] we have a^^\v) — Ca^^\v) then 

that is, the parametrization is trivially in the sense that it can be made by scaling the 
variable by a multiplicative factor in (one of) the two patches. If we are in such a case, as for 
parametrizations, one can easily conclude 5r([0,1]) x 5r([0,1]) = Vp. If we are not in such 
a special configuration, then Theorem states that we can independently approximate the trace 
and transversal derivative of a function by isogeometric functions, but for /i-refinement the 
approximation error does not converge to zero. 


5.2. General two-pateh geometry 

In this section we study the approximation property of V^, for two-patch geometry parametriza¬ 
tions beyond analysis-suitable -continuity. We consider a specific case where is the identity 
mapping while G [5r and select 


Q/(^)(p) = 1 

a^^\v) = DuF^^\0,v)-ei G ^^([O,!]), 

p(^\v) = 0 

= DuF^^\0,v)-e2 e 5 ?([ 0 , 1 ]). 


Then d = ei and 

= vf = 


Equation (29) simplifies to 



such that 0 G V 



0 


d(^ 
’ dx 




oFo = [go, 9 i] 

g^^\u, v) = 5fo(^^) + gi{v) u + O(m^), 

g(^\u,v) = go{v) + {(3^^\v)gQ{v)+ a^^\v)gi{v))u + 0{u^). 


(33) 


(34) 


(35) 


The following statement gives a full characterization of the space Vp. We use ^([0,1]) to indicate 
the null space 0, and we recall that = 0 is not allowed. 
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Theorem 3. Let Vi = U he the identity and 

Assume G 5r'^([0,1])^ with 0 5r'^~^([0,1]). 

If = 0; then 


a(^\ /3(^), /3(^) as in ([ 33 I). 


Vi=5?([0,l])x 5^(10,1]). (36) 

If instead G 5r^([0,1]) with 0 Sr^ ^([0,1])^ then 

([0,1]) X 5r'^-([0,1]). (37) 

Proof. Since is the identity we have Vp = Vp C 5^ ([0,1]) x 5r ([0,1]). 

Consider first the case = 0. It is easy to see that 


ViD5,^([0,l])x5r^-([0,l]). 


Indeed, given any [go^gi] G 5r([0,1]) x Sr ^"^([O,!]) we can find g^^\u^v) G such that 

(35) holds. In particular, take 


g^^\u,v) ^ g(){v) P a^^\v)gi{v)u. 


{R)( 


Moreover, 5?([0,1]) X Sr ^"^([0,1]) is equal to the space Vp due to the second condition in ( [3^ , 
which is 

g^^\u,v) = go{v) + a^^\v)gi{v) u + 0{u^) 

for = 0, since g^^\u, v) G <Sr forbids gi to be of a polynomial degree higher than p — 

The second case, G <Sr^([0,1]) with ^ 5f^ ^([0,1]), is similar. Again, we have Vp D 
^mn{p,p P/3+^i^QQ^ x 5?([0,1]). Indeed, given any 

bo,5i] e ([0,1]) X ^^([O,!]) 

we can find g^^\u^v) G Sripf^^^) such that ( [35| ) holds. This time, take 

g^^\u,v) = go{v) + {P^^^v)go{v) + a^^\v)gi{v)) u. 

As before ^^^^^([0, Ij) x <Sr~^“([0,1]) is eqnal to the space Vp due to the second condition 

in (35), that is 

3 g^^\u,v) = 5fo(^^) + {l3^^\v)go{v) + a^^\v)gi{v)) u + 0{u^). 


This concludes the proof. 


□ 


Remark 7. In the setting of this seetion, F^^) being the identity and F^^) G Sr(Li^^^) x Sr(Li^^^); 
we always have Pa ^ p, Pg ^ p, whieh follows from (33). Then Theorem^ guarantees 

Vf = VI D 5 ^( 0 , Ij) X 5 ° i([0, Ij) ipo. 


This is in aeeordanee with the well known faet that isoparametrie funetions reproduee all linear 
polynomials. 
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We remark that any restriction on Vp as stated in Theorem affects the approximation prop¬ 
erties of V^I^(L) C and consequently of itself. For the sake of convenience, this is 

summarized in the corollary below. 

Corollary 1. Optimal order of convergence for h-refinement can not be achieved if deg > 1 
or deg/?^^) >1. In particular locking, i.e. no convergence, occurs for degcr^^^ > p — r or 
deg > p — r. 

Some examples will be considered and further analyzed in Section 

6. isogeometric spaces on surfaces 

The result of the previous sections can be extended to a a multi-patch surface 


We can set up an isogeometric function space over the surface 

V =!</>: ^ R such that 0 o € SP{Q),i = 1,..., At} , (38) 

and again = V fl and = V H As for the planar case, : Q though 

is now a surface patch. More important, for the function space to be well defined, the 

surface itself needs to be i.e., the surface needs to have a well defined tangent plane in every 
point. For simplicity, we focus on a single interface, that is a two-patch geometry, where each patch 
is parameterized via 

F(^) : [-1,0] X [0,1] = n^^'> C 

F^^) : [0, 1] X [0, 1] = ^ C 

with G (5r for S G {L,R}. We ask the parameterization to be i.e., there exist 

: [0,1] ^ M, : [0,1] ^ M and /3 : [0,1] ^ M such that \/v G [0,1], a^^\v)a^^\v) > 0 and 

a^^\v)DuF^^\0,v) - a^^\v)DuF^^\0,v) + f3{v)D^Fo{v) = 0. (40) 

Then, the surface gradient of a function 0 : ^ M can be computed as follows. First, the surface 

is extended to C by defining 

x]-e,e[^J4f) (41) 

with 

{u, re) = {u, v) + reN^^) {u, v ), (42) 

where is the unit normal vector of Now, we define the extension $ of the function 

(f via for S G The surface gradient of (j) is then given as the 

three-dimensional gradient of the extension $ in restricted to the surface fl, i.e. 

Vq 0 = v$|n, (43) 
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where 


(44) 


v$ O on 


By construction, the gradient is tangential to the surface. 

Proposition can be generalized to surface domains. For the sake of simplicity we consider a 
simplified case, by assuming that can be projected onto the (x,y)-plane without self intersec¬ 
tions. Note that this is not a limitation of the concept but facilitates the following propositions. 
Then we have 

f(^) ^ op(^))^ (45) 

where P^*^) = is the planar parameterization of the projected surface and fs is an 

(isogeometric) function from the planar projection = p('^)(ri('^)) to M. Then the following 
result is a direct corollary of Proposition see also [I6] . 


Proposition 4. An isogeometric function 0 G V as in (38) belongs to if and only if its graph 
as a 2-dimensional surface in is -continuous at each interface H 


The definition for a G^-continuous graph of an isogeometric function (f is formally the same as 
in Definition!^ however now the graph of 0 : D ^ M is a two-dimensional surface in M^. As for 
planar domains considered in Section the coefficients for the condition are still determined 
by the parameterization in fact by two of its three components. In our simplified case, we are 
in the following situation. 


Proposition 5. Consider surface patches 


and functions g^^\ g^^^ fulfilling Definition^ Then the coefficients 
given by a^^\v) = , for S G and P{v) = ^{v)P{v), where 


and I3 in (40) are 


a^^\v) = det 


DuP^^\0,v) D^p(^^0,v) , 


(46) 


^{v) = det 


DuF(^\0,v) DuF(^\0,v) , 


(47) 


and 7 : [0,1] —> M is any scalar function. In addition, '-f{v) ^ 0 if and only if a^^\v)a^^\v) > 0. 
Moreover, there exist funetions I3^^\v), for S € {L,R}, sueh that 


P{v) = - a^^\v)F\v). 


(48) 


Again, we omit the details of the proof, as it is a direct consequence of Proposition 
In analogy to the planar case, to guarantee optimal approximation properties of the isogeo¬ 
metric space we ask some structure for the trace and derivative trace of isogeometric functions, 
by restricting the geometry parametrization to the class of analysis-suitable G^-continuous func¬ 
tions. The definition is omitted being exactly the same as Definition In the case of ( [4^ , the 
parametrization P^*^) of the planar projection D of the geometry D needs to be AS G^. Furthermore, 
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the third component of as well as the function have to fulfill the same conditions, for 
S G This means that Remarks |3]-[^ extend to surfaces: isogeometric spaces over and 

AS multi-patch parametrizations fit in the isoparametric framework. Eventually, (j) G if 

and only if 0 G V^(ri), where 0 o =00 Then, results of Section apply directly to AS 


G^ surfaces, at least when the geometry parametrization fulfills (45). 


Remark 8. As a consequence of Proposition\^ G^ isogeometric spaces over G^ multi-patch surface 
parametrizations fit in the isoparametric framework. 


7. NURBS spaces 

The results presented in Sections |3]-[^ are not limited to spline spaces, but, for example, are 
generalizable to NURBS. One possibility is to generate rational planar geometries from AS G^ 
surface patches. A rational patch 

(49) 


F® : h 54® C 


with 




can be interpreted as a surface patch in M^, with 


(50) 


(51) 


when transformed to homogeneous coordinates. Recall, that two points F, F' in homogeneous 
coordinates correspond to the same point in Euclidean space, if there exists a A ^ 0 such that 
F = AF'. An AS G^ surface with surface patches (51) generates a rational AS G^ geometry via 
(50). Therefore this construction fits into the framework of Section]^ For the sake of brevity, we 
only present one example (Figure]^ in the next Section 8.2.2 and postpone further studies. We 
want to point out that not all multi-patch NURBS geometries fit into this framework. The surface 


representation in (51) may not be G^ (or not even G^) but the rational patches representing the 
graph still join G^. Such a configuration is given for the circle presented in Figure [T^ However, in 
that example the graph surface is not analysis-suitable G^. 


8. Numerical tests 

In this section, we present numerical tests which illustrate the previous theoretical results. The 
simulations have been obtained by the isogeometric analysis library IGATOOLS, see m for an 
introduction. For related numerical studies see [23 EH Ea ED]. 
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8.1. Model problem 

We consider the following bilaplacian problem on C where the unknown is denoted by w 
and the data by /, 

A‘^w = f 
< w = 0 

Vw •n = 0 

In our tests the data / is selected in order to have an analytic exact solution w. Here n is the unit 
normal vector to the boundary d^. Let 


n, 

on, (52) 

on. 


Uo = {v ^ = 0 on dQ and Vv • n = 0 on 512}, 


the variational form of (52) is, find w G ZYq, such that 

/ AwAvdxdy= / fv dxdy^ 

Jn Jn 


(53) 


In what follows, we consider two types of geometries, the ones which are analysis-suitable 
and the ones which are not. Let us start with the analysis-suitable geometries. 


8.2. Analysis-suitable G^ geometries 

We consider five different analysis-suitable G^ geometry mappings with different numbers of 
patches. For all cases, we select the degrees p = 3,4, 5 and regularities 1 < r < p — 1. 


8.2.1. Two-pateh geometry (L-shape) 



Figure 4: Two-patch L-shaped domain (left), a solution affected by locking for p = 3, r = 2 (center), and a correct 
numerical solution for p = 3, r = 1 (right). 

We start with the L-shaped geometry consisting of two patches. The L-shaped domain is given in 
Figure]^ (left), where the red edge is the common edge of both patches. Here the parametrization 
of both patches is bilinear. This is obviously an AS G^ geometry (recall Proposition]^. Using 
Theorem [T| the optimal convergence is achieved if and only if r < p — 2, in agreement with the 
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results of Section If we focus on degree p = 3^ locking is evident when the regularity equals 
to r = ^ — 1 = 2, see Figure]^ (center). In particular, the solution on the interface line equals to 0. 
In order to circumvent locking, we decrease the regularity, see Figure(right). Figure]^ gives 
the convergence curves for degrees p = 3,4 and 5. We obtain expected results for all degrees, i.e. 
we have convergence (which is optimal) if and only if r + 2 <p. 


J2 gJ-J-QJ-g 



Level of refinements 


Reference 

- - - p = 3 and r = 2 

p — 3 and r = 1 
Reference 

-, -1 - p = 4 and r = 3 

- - - p = 4: and r = 2 

p = 4 and r = 1 
Reference 

... p = 5 and r = 4 
... I - p = 5 and r = 3 
... p — 3 and r = 2 
p = 5 and r = 1 


Figure 5: Convergence results for the L-shaped domain. 


8.2.2. Multi-patch geometries (triangle, quarter of a circle and rectangle) 



Figure 6: Three-patch triangle (left), a solution affected by locking for p = 3, r = 2 (center), and a correct 
numerical solution for p = 3, r = 1 (right). 

In what follows, we consider cases with three or more patches, starting with an example with 
three bilinear patches forming a triangle (see Figure]^ (l^ff)) another with three bi-quadratic 
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Figure 7: Three-patch quarter of a circle (left), a solution affected by locking for p = 3, r = 2 (center), and a 
correct numerical solution for p = 3, r = 1 (right). 

patches forming the quarter of a circle (see Figure (left)). Both are AS geometries. Note 
that the quarter of the circle is composed of NURBS patches, and is obtained following Section 
construction, from a geometry parametrization that is an AS surface in homogeneous coordi¬ 
nates. The details of this construction are presented in the appendix. Figures (center) and 
(center) show G^ locking which appears with p = 3 and r — 2 that we can circumvent by reducing 
the regularity, see Figures]^ (right) and (right). The expected convergence orders are given in 
Figure 


L2 


errors 


L2 


errors 



Reference 
p — 3 and r = 2 
p = 3 and r = 1 
Reference 
p = 4: and r = 3 
p = 4 and r = 2 
p — 4 and r — 1 
Reference /i® 
p — b and r — 4 
p — b and r — 3 
p = b and r = 2 
p — b and r — 1 


1 2 
Level of refinements 




Reference 



- p 

— 3 and r 

= 

2 

— p 

^ 3 and r 

- 

1 


Reference 



.,.p 

— 4 and r 

= 

3 

- p 

— 4 and r 


2 

— p 

— 4 and r 

- 

1 


Reference 



- p 

— 5 and r 

= 

4 


— 5 and r 

- 

3 

- p 

^ 5 and r 


2 

~p 

— 5 and r 

- 

1 


Figure 8: Convergence results for the triangle (left) and the quarter of circle (right). 

Next we consider a rectangle composed of four patches, see Figure (left). As shown in this 
figure, we consider a particular case where two interfaces are collinear. This configuration is among 
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Figure 9: Four-patch rectangle (left), a solution affected by locking for p = 3, r = 2 (center), and a correct 
numerical solution for p = 3, r = 1 (right). 


the ones analyzed in [HI Section 11.2.3]. The results are presented in Figure for degree p = 3 and 
regularity r = 2 (center) and r — 1 (right). Figure [TT| (left) gives the convergence orders, which are 
as expected. 



Figure 10: Five-patch simply-connected domain with smooth boundary (left), a solution affected by locking for 
p = 3, r = 2 (center), and a correct numerical solution for p = 3, r = 1 (right). 

The final and most relevant example of AS geometry is reported in Figure [To| (left). This is a 
five patch decomposition of a simply-connected domain with smooth boundary. The parametriza- 
tion of each patch is bi-quadratic, and the domain boundary is C^. Given the boundary control 
points and parametrization, the interior contol points have been selected in order to fulfil the AS 
conditions. Unlike [32], here not only the mesh curves, i.e., patch interfaces, but also their 
parametrization need to be chosen properly. We refer to the Appendix for the complete descrip¬ 
tion. The results are presented in Figure for degree p = 3 and regularity r = 2 (center) and 
r = 1 (right). Figure [TT| (right) gives the convergence orders, which are again as expected. 

8.3. Non-analysis-suitable G^ geometries 

In this section we study two different examples of non-analysis-suitable G^ geometries. First we 
consider a two-patch parametrization of a rectangle, where we compare a distorted parametrization 
with the undistorted identity mapping, see Figure As a second example we consider a five-patch 
circle. 
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... p = 3 and r = 2 
p = 3 and r = 1 
Reference 

.,., - p ^ 4 and r = 3 
... p — 4: and r — 2 
p — 4 and r — 1 
Reference /i® 

. I . p — 3 and r — 4 
.,., - p ^ 5 and r = 3 
... p = 3 and r = 2 
p — b and r — 1 



Reference 

... p ^ 3 and r = 2 
^ 3 and r = 1 
Reference 

... p = 4 and r = 3 
. I., - p — 4 and r = 2 
p — 4 and r — 1 
Reference 

... p — 5 and r = 4 
.,., - p ^ 5 and r = 3 
p = b and r = 2 


Figure 11: Convergence results for the four-patch rectangle (left) and the five-patch simply-connected domain with 
smooth boundary (right). 


8.3.1. Two-patch geometry (quadratically distorted rectangle) 



Figure 12: Non-AS rectangle (left), two solutions affected by locking for p = 3, r = 2 (center-left) and p = 3, 

r = 1 (center-right), as well as a correct numerical solution for p = 4, r = 1 (right). 


We consider a rectangle domain [—1,1] x [0,1] formed by two patches, where the left one is 
linear (identity) and right one is quadratic in the horizontal direction and linear in the vertical, 
see Figure [T^ (left). This case illustrates the configuration covered by Theorem]^ and Corollary 
As shown in Figure [T^ with degree p = 3, (7^ locking is manifested for both r — 2 and r = 1, see 
respectively the second and the third columns of the figure. In order to have convergence, a degree 
of at least p — A has to be selected, see Figure [T^ (right). 

However, as anticipated in Corollary [Tj we cannot expect optimal convergence. This is a direct 
implication of Theorem]^ which for this case (pg, = 2 and (3^^^ — 0) states that for all 0 G 


d(j) 

dx 


G5r^-(F) =5f(F). 


(54) 


Assume that the exact solution w is smooth enough and let be the numerical solution. Then, us¬ 
ing (54), together with the usual approximation estimates in Sobolev norms and the trace inequality 
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for Sobolev spaces, gives 


'^approx' — 


dw dwh 


^ ^ CtraceW"^ '^h\\H^[Q)') 

H2(r) 


dx dx 

instead of the optimal order of convergence, that is when measuring the error in the i7^(f])-norm. 



Figure 13: Two-patch rectangle with identity mapping (left), two-patch rectangle with a quadratic distorted patch 
(right). 


In Figure we report the convergence results and, due to the specificity of this case, both 
L? and norms are plotted. While locking is easily recognized, from these numerical tests 
it is difficult to measure the expected sub-optimality of the asymptotic behavior. This is likely 
a numerical artifact due to the imposition of the constraint in our implementation, which is 
discussed in Section 8.4 We further analyze this example in Figure 15, where we compute, for 
p = 4 and r = 1, the error only on the left patch, that is [—1,0] x [0,1], and compare the results for 
this geometry and for a reference geometry formed by the identity mapping on both patches (see 
Figure . Note that on the left patch both parametrizations are the identity mapping. 


^2 gJ-J-QJ-g 




J_^^^_L 

0 12 3 4 


Level of refinements 


Figure 14: Convergence results for the non-AS two-patch rectangle, with error in (left) and (right). 
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Level of refinements 


Level of refinements 


Figure 15: Comparison of convergence between the two geometries given in Figure for p = 4 and r = 1; the error 
is computed only in the left subdomain [—1,0] x [0,1] in (left plot) and (right plot). 



Figure 16: Five-patch circle (left), two solutions affected by locking for p = 3, r = 2 (center-left) and for p = 3, 
r — 1 (center-right), as well as a correct numerical solution for p = 4, r = 1 (right). 


8.3.2. Multi-patch geometry (circle) 

In the last example, we study an exact circle composed of five patches, given in Figure 
(left). Here, we are interested in testing rational parametrizations that are beyond the framework 
presented in Sectionj^ since their homogeneous representation is not an analysis-suitable surface. 
Even though our theory does not apply, the numerical results obtained are consistent with our 
findings. For degree p = 3 the numerical solution suffers of locking, as one can see in Figure [T^ 
for regularity r — 2 (center-left) and r — 1 (center-right). For degree p — A and regularity r = 1, 
see Figure [T^ (right), we observe convergence to a solution. Figure [T7| (right) gives the convergence 
behavior for degrees p = 3,4 and 5. Sub-optimality in all situations is manifested. 
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Figure 17: Convergence results for the five-patch circle. 


5 . 4 . Numerical implementation of continuity 


In this section, we describe the numerical implementation of continuity that we have used 
in order to obtain all the numerical examples given previously. Let A and 6, be the matrix and the 
right and side of the system corresponding to the variational formulation (53), where no boundary 
condition or continuity condition are included. Furthermore, let Ci be the constraint matrix in 
symmetric form, i.e. 


where n is a normal unitary vector and [-J is the jump at the patch interface F. Let Nq be the change 
of basis matrix from the fully unconstrained space to the subspace fulfilling boundary condition and 
continuity. Nq can be seen as obtained from the identity I Nunc ~ where Nunc is the number 
of degrees of freedom without constraints - by removing the columns with index corresponding 
to the degrees of freedom of the boundary conditions and by summing the columns of degrees of 
freedom that are shared on F. Since it is not trivial to compute a continuous basis analytically, 
we operate numerically by computing the null-space 


Ni =null(N^CiNo). 


Then we solve the following problem: find 2 ; such that 


NfN^ANoNi z = NfN^6, 


(55) 


(56) 


and obtain the solution in the unconstrained initial spline basis as x = NqNiz. 
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The numerical computation of (55) is a hard task for non AS geometries. Indeed, the 
non-zero eigenvalues of N^CiNq are not well separated from the eigenvalues that are numerically 
zero (close to machine precision), as shown in Figure 18 for the distorted rectangular domain in 
Figure 12 (left), with p = 4, r = 1. For the next refinement step, the distinction may not be 
possible anymore. By comparison. Figure shows the eigenvalues for the L-shaped domain with 
the same degree and regularity. In this case, it is easier to distinguish between numerical zeros 
(around 10“^^) and non-zero eigenvalues (above 10“^). The two configurations are structurally 
equivalent, i.e. the topology of the control point grid is the same. 


Eigenvalues 



Figure 18: Eigenvalues of the domain given in Figure(left) for p = 4 and r = 1, and different /i-refinement levels. 


Eigenvalues 



•ref = 4 
•ref = 3 
•ref = 2 
•ref = 1 
ref = 0 


Figure 19: Eigenvalues of the L-shaped domain given in Figure]^ (left) for p = 4 and r = 1, and different /i-refinement 
levels. 


9. Conclusions 

In this paper, we have studied C^-smooth isogeometric function spaces over multi-patch geome¬ 
tries. We have considered geometry parametrizations composed of multiple patches that meet at 
patch interfaces. As it is common for isogeometric methods, the geometry parametrization is given 
initially and considered fixed. From that, the geometric continuity conditions are derived. The 
same conditions need to be satisfied for the graph parametrizations of the isogeometric functions 
in order to obtain a C^-smooth function space over the given geometry. 

We have studied how these conditions affect the traces and the transversal derivative of the 
isogeometric function spaces along the patch interfaces.Indeed, if the trace or transversal derivative 
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at the interface is over-constrained, the approximation order of the isogeometric function space is 
reduced. We identify a class of configurations that allows for optimal approximation properties 
of isogeometric spaces, so called analysis-suitable (AS G^) geometry parametrizations. We 
show numerically that AS G^ geometries indeed allow for optimal approximation. We have not 
developed a complete approximation theory with classical error estimates, but this will be the topic 
of a future paper. Parametrizations that are not AS G^ cause suboptimal order of approximation. 
In the worst case the convergence under /i-refinement is prohibited, a behaviour that we have named 
G^ locking. 

We have addressed mainly the case of planar B-spline geometries, but we have briefly discussed 
the generalization to surfaces and NURBS patches. All the results are supported and confirmed by 
numerical tests for various degrees and orders of regularity of the spline space. We have numerically 
solved a bilaplacian problem over several AS and non AS G^ geometries, by a standard Galerkin 
approach. As we pointed out, the numerical implementation of the G^ conditions poses some 
non-negligible difficulties for complex configurations of non AS G^ geometries. 

An important question that remains to be studied in more detail is the flexibility of analysis- 
suitable G^ geometries. As we have shown, the AS G^ class contains bilinear patches but extends 
to more general configurations. We formulate the problem of flexibility in the following way: Given 
a collection of boundary curves, is it possible to find patches that interpolate the boundary curves 
and that form an AS G^ geometry? For piecewise linear boundary curves, the problem can be 
solved by using bilinear patches. In Figure]^ we have shown the interesting example of a piecewise 
biquadratic AS G^ parametrization of a G^ simply-connected domain. The extension to arbitrary 
degrees and topology deserves further investigation. 

The construction for AS G^ surfaces could be more difficult. In the planar case, the interior 
parametrization of the patches may be modified in order to achieve analysis-suitability. This is not 
feasible for surfaces, as the surface itself changes if the parametrization of the interior is changed. 
AS G^ constructions are possible if the surface is given as the image of a planar domain. For 
the general setting, an explicit construction remains an open problem to be considered in future 
research. 
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Appendix A: AS parametrization of a quarter of circle 


In this appendix we discuss in detail the construction of the geometry parametrization of the 
example shown in Figure in Section [8. 2. 2[ This is based on the NURBS setting of Section and 
on the ideas in [8] as well as m Section 3.4]. Each patch of the three-patch geometry is constructed 
from a combination of a bilinear mapping 

B« : [0,1]2 ^ c A = {(n, v) G [0,1]2 : u + v< 1} 
onto a quadrilateral within a reference triangle A and a global mapping 

G : A ^ ^ C 


from the reference triangle to the surface patch Q in homogeneous coordinates. These mappings 
and corresponding domains are shown in Figure 


[ 0 , 1]2 


B« 



Figure 20: Mappings and G and corresponding domains. 


Here the mapping G is a triangular Bezier patch of total degree p = 2 with 

_ ^ j\ h\ 

G{s,t)= ^ s" 0^-s-tf 

i+j+k=2 

with control points 

50 , 0,2 = ( 0 , 0 , 1 )^ 50 , 1,1 = ( 0 , 2 \/ 2 )^ 50 , 2,0 = ( 0 , 1 ; 1 )^ 

51,0,1 = (V2,0, 2V2f 51,1,0 = (2^2, 2V2, 2V2f 
52,0,0 = (1,0,1)'^ 

in homogeneous coordinates. In Figurethe dashed blue lines represent the corresponding control 
point grids in Cartesian coordinates. Each quadrilateral and is formed by one corner 

point of the triangle A, two adjacent edge midpoints as well as the center of gravity (1/3,1/3)^. 
The bilinear mapping is (up to rotations of the parameter domain) completely determined by 
its image By construction, the mapping = G o is a rational bi-quadratic function in 
each component. 

Using this construction, the theory developed in Sections and can be applied to the example 
in Figure 
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Appendix B: AS parametrization of a smooth simply-connected domain 


In the following we give a description of the domain depicted in Figure The domain is 
composed of five bi-quadratic patches and the boundary is C^-smooth. The central patch is just 
the shifted unit square parameterized by 

F^'^\u,v) = ^ 

The top patch is parameterized by 

y2v^u-u^) + l(l + ^v-^v^) ) 

and for the left, bottom and right patches are given by 
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